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Abstract 

We found that spin-one new light particle exchanges are strongly 
bounded by high-energy and small momentum transfer np elastic scat- 
tering data; the analogous bound for a scalar particle is considerably 
weaker, while for a pseudoscalar particle no bounds can be set. These 
bounds are compared with the bounds extracted from low-energy 
n — Pb scattering experiments and from the bounds of ir° and K + 
meson decays. 



1 Introduction 

The Standard Model of three fundamental forces describes interactions of 
elementary particles very well. While the electromagnetic force has a long 
interaction range, the short radius of the "weak force" (~ 1/1000 fermi) 
is determined by the heavy masses of mediating W and Z bosons (~ 100 
GeV). The QCD forces are typically contained within a confinement radius 
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of ~ 1 fermi. The effects of the long-range gravitational force can usually be 
neglected in elementary particle scattering experiments. 

The search for new forces of nature is a major goal of experiments at 
high-energy colliders. Rare transitions and decays of fundamental particles 
can also shed light on new interactions. These experiments are probing for 
new forces at distances shorter than 1/1000 fermi. However, the existence of 
new forces at distances larger than the confinement radius of a nucleon (~ 1 
fermi) can also be probed and constrained by sensitive experiments. 

Experimental searches for and limits on these new forces of nature can 
be pursued in two directions: (a) as deviations from the Newtonian law of 
gravity where the new force is expressed as a modification of the 1/r 2 law, 
usually by an additional Yukawa term that can be parameterized with two 
parameters a, the relative strength of new interaction, and A, the character- 
istic radius of the interaction. Applied in the analysis of experimental data 
at macroscopic distances down to ~ a micrometer, this ansatz describes the 
possible deviations from classical gravity; (b) as a quantum field theory de- 
scription of the interactions (excluding gravity) in a covariant form, which 
can be expressed in the lowest perturbation order through the coupling con- 
stant g and the mass of the exchanged particle mediating the interaction /i. 
Covariant forms that can be consistently considered in this description are 
scalar (S), pseudoscalar (P), vector (V), and axial vector (A). We can argue 
that higher spins of the intermediate particle should not be considered since 
they lead to non-renormalizable theory. The particles that mediate the new 
force could be absent from the spectrum of known particles [1] due to their 
small mass and coupling constant or due to some other reason that is helping 
them avoid detection. In any case, if these particles are not observed, direct 
experimental limits on their existence in terms of g and \x are required. 

In Section 2 of the present paper, we reanalyze the experimental small- 
angle np-elastic scattering data at high energy [2] in terms of bounds on the 
existence of new forces expressed as S, P, V, or A covariant interactions. In 
Section 3, we examine bounds that can be obtained from lower energy data. 

2 Bounds from high-energy np scattering 

The data for small-angle np-elastic scattering at high-energy were obtained 
in the NA-6 experiment [2] performed at CERN SPS a quarter century ago. 
Incident neutron energy in the experiment was 100-400 GeV, while the square 
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of the 4-momentum transfer \t\ was varied in the range 6-10 -6 to 5-10 -1 GeV 2 . 
The data of this experiment are consistent with extrapolation of the hadronic 
amplitude from higher \t\ values, while at \t\ < 10~ 4 GeV 2 the differential 
cross-section rises due to Schwinger scattering, which is the interaction of 
the neutron's magnetic moment with the Coulomb field of the proton or 
electron. The purpose of NA-6 [2] was to measure hadronic interactions at 
high s in the region of momentum transfer (~ \t\ < 10 -2 GeV 2 ) that was 
usually inaccessible in the scattering of charged hadrons due to Coulomb 
interactions. This is the region where the effect of a new force mediated by 
a light particle may be present. 

Figure [l] (similar to Fig. 16 from [2]) demonstrates that np elastic scat- 
tering data in this experiment are well described by the following formula: 

dcr . ,, , „ / ak n \ 2 n . 
- = ,4exp[ W ] - 2 (— j - , (1) 

where A = (79.78 ± 0.26)m6/GeV 2 and b = (11.63 ± 0.08) GeV~ 2 were 
determined from the fit to the data (data are taken from Table 7 in [2]), m n is 
the neutron mass, and k n = —1.91 is the neutron magnetic moment in nuclear 
magnetons. The factor of 2 in the Schwinger term, as will be discussed later, 
accounts for the scattering of the neutron's magnetic moment on the proton 
plus an incoherent contribution of scattering on electrons (gaseous hydrogen 
was used in [2] as a target). Smaller effects due to neutron polarizability 
are not included in the description of the data. This description (JTJ works 
rather satisfactorily with x 2 =41.5 for 31 degrees of freedom. We will refer 
to this description as the "zero model" since no new force contributions are 
included here. 

Although Quantum Chromodynamics does not provide a detailed theo- 
retical description of the hadronic elastic scattering at small \t\, i.e. at large 
impact parameters, hadronic scattering has been studied experimentally in 
great detail in the past and was phenomenologically well understood, e.g. in 
the framework of Regge models. The description of elastic data by a single 
exponent was a general universal feature of hadronic scattering observed at 
low \t\ in the region where it was not obstructed by Coulomb scattering (for 
example, see [3] and also the comparison with other experiments in [2]). This 
justifies, at a phenomenological level, our choice of the hadronic scattering 
description with a single exponent. However, in an attempt to improve the 
description of the data [2] , we have tried several alternative modifications of 
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Figure 1: Elastic differential neutron- proton cross-sections measured in experi- 
ment [2]. For comparison, the \t\ region measurable in pp scattering is shown with 
the effect of the Coulomb interaction indicated by the dashed line. 
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the exponential term in the "zero model" Q involving additional parameters, 
including a quadratic term in the exponent and the sum of two exponents. 
In all of these cases, \ 2 P er degrees of freedom was slightly increased demon- 
strating that more complicated modifications of the "zero model" are not 
statistically justifiable. 

We describe the contribution of a new interaction in the following way: 
Let us suppose that a new light particle with mass /i exists which interacts 
with the neutron and proton with couplings g n and g p correspondingly. As- 
suming scalar, pseudoscalar, vector, and axial vector couplings of this particle 
with nucleons, we obtain the following addition to expression (jl|: 

dai |A1 2 • FF 

dt [9 ^ )lnew ~ lQns(s-Am^) ' 1 ] 

where s = (p n + p p ) 2 is the invariant energy square and m is the nucleon 
mass. We parameterize the hadronic form factor's contribution as: 

FF= (l-t/A2)8 ' (3) 

which comes from a 1/g 4 decrease of the nucleon form factor, and we set A 
equal to the mass of the lightest meson resonance with appropriate quan- 
tum numbers (7/ in the case of pseudoscalar). Finally, we use the following 
amplitude squares for different couplings: 

\A S \ 2 = 7 ^K^{^-tf , (4) 



(t - ' 

(t-/i 2 ^ 2 



\M 2 = , (5) 



|A y | 2 = -i^-[ S 2 -4m 2 S + 4m 4 + S t+^ 2 ] , (6) 

A A 2 = 7, y \ 2 [s 2 + ±m 2 s + 4m 4 + st + -t 2 + + — , 7 

where coupling constants gf = g p g l n - 

It is quite natural to suppose that a new light particle's couplings with 
nucleons originates from its couplings with quarks. In this case, (6) and 
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(7) are modified. For the vector exchange, the induced magnetic moment's 
interaction term should be added to the scattering amplitude. Since its 
numerator contains momentum transfer divided by m^v, which in considered 
kinematics gives a factor much smaller than 1, we can safely neglect it and 
use (6) in what follows. The case of the axial vector exchange is more delicate 
and discussed in detail in the Appendix. 

We can now turn to the discussion of other features of the np elastic 
scattering amplitude. Though the strong interaction amplitude cannot be 
determined theoretically from the first principles, our confidence that 1/t 
dependence is absent in strong interactions for \t\ < opens the road to 
bounding the light particle exchange if its mass is smaller than that of the 
7r-meson. Experimental data at \t\ < m\ matter for our bounds, which makes 
the precise value of A in the expression for FF not important, since in the 
relevant domain of \t\ the form factor is close to 1. For the same reason, no 
form factor is introduced for the Schwinger term in ([I]). 

For each fixed set of parameters gf and /x describing the possible contri- 
bution of a "new force" , we are fitting the experimental distribution with a 
combined function (JXj) — |— (|2j) , where parameters A and b describing the stan- 
dard hadronic contribution are free. Then the maps of A, b, and the minimum 
values of x 2 are composed as functions of gf and fi. Analyzing the x 2 map, 
we determined the level of x 2 BO above which parameters of the "new force" 
become incompatible with experimental data at a confidence level (C.L.) 
greater than 90%. At this level we also examined and ensured that parame- 
ters A and b remain within the 90% C.L. close to those in the "zero model". 
In this way, we can ensure that the "new force" contribution does not sub- 
stitute for the standard hadronic plus electromagnetic contributions in the 
description of the data. 

Figure [2] shows, for comparison, fits to the data for the "zero model" 
and for several excluded models for the new vector particle exchange with 
parameters slightly beyond the excluded limits for /zy and g\. 

Two comments need to be made on formulas (|4])-([7]): (a) the amplitudes 
with the exchange in the t-channel of a point like particle with spin a de- 
pend on s as s a . That results in an amplitude behavior of s° for scalar and 
pseudoscalar and of s for vector and axial vector particles. This property 
of high-energy scattering amplitudes would allow us to determine the value 
of the spin of the "new physics" mediator; (b) the pseudoscalar exchange 
vanishes at t — 0. These comments explain why we will get the strongest 
bounds on gA and gv, a weaker bound on g$, and no bound on gp. 
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Figure 2: Several fits to the experimental data of [2j: long-dash line - single expo- 
nent without the Schwinger contribution; solid line - the "zero model" description 
of Q; dotted line - the "zero model" plus the new vector particle contribution 
with fx = 1 MeV and g 2 = 0.0015; short-dash line - the "zero model" plus the new 
vector particle contribution with \i = 10 MeV and g 2 = 0.005; dot-dashed line - 
the "zero model" plus the new vector particle contribution with fi = 40 MeV and 
g 2 = 0.025. 



Before presenting fit results, let us explain why we neglect the interference 
of a new particle exchange amplitude with the strong amplitude and with 
the photon exchange (Schwinger) amplitude. The strong amplitude is almost 
entirely imaginary in the energy domain studied in j2] (\Re/Im\ < 0.1 [1]). 
That is why it does not interfere with the real amplitude of a point-like new 
particle exchange. Interference of the strong amplitude (as well as the new 
force amplitude) with the Schwinger term is negligible since the interference 
term is constant at t = Q unlike the square of the Schwinger amplitude, 
which contributes significantly at small \t\ because of 1/t behavior. 

The mass of a light particle \i was bounded in our fits to be below 100 
MeV, and the range of the coupling constants varied depending on the par- 
ticular model. Compilation of the bounds obtained from the y 2 limit for P, 
S, V, and A models in coordinates g 2 versus /x is presented in Figure [3} 

In the next step of analysis, we checked that for each model, fitted pa- 
rameters A and b corresponding to the boundary of the excluded domain 
of g 2 and /i must not deviate from their "zero model" values by more than 
1.28cr, where a is the corresponding error of parameters A and b from the 
"zero model." We found that these conditions are satisfied if \x < AOMeV for 
the S, V, and A-models, and cannot be satisfied for any value of /i for the 
P-model. Thus, our bounds on the coupling strength g 2 shown in Figure k| 
should only be referred to in these validated domains (indicated in Figure w 
by brackets). No consistent limit can be set for the P-model. In addition, we 
should notice that the very high value of gp obtained from the x 2 analysis 
for the P-model makes our perturbative approach of formula ^ not valid. 
We should conclude therefore that the experimental data [2J do not provide 
any limit for the pseudoscalar exchange. 

The factor of 2 in the Schwinger term of the "zero model" in ([I]) is coming 
from both n — p and n — e scattering and is an estimate of equal contribution 
from both. However, n — e scattering occurs in a different kinematical range 
and the event selection criteria in [2] could suppress the detection of electrons. 
Consequently, we varied the factor in the Schwinger term of the "zero model" 
([I]) from 1 (no n — e contribution) to 3 (double n — e contribution) and found 
in analysis that this variation was not very significant, changing our limiting 
value for g 2 by ± 8% for a fixed value of \i. 

One can notice that since the average s ~ 540 GeV 2 in experiment [2], 

lr The Schwinger part of the interference term contains q^/q 2 , q 2 = t, which multiplies 
= (j>i — P2)fi cancelling the l/q 2 enhancement, or (pi + P2)^ giving zero. 



8 




Figure 3: Compilation of the upper bounds obtained in the current analysis in 
terms of g 2 and fi at a 90% C.L. No limit can be set for the pseudoscalar exchange. 
Brackets indicate the interval of mass /x where the analysis was validated (see text). 
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vector (Ay) and axial vector (A a) amplitudes, as follows from Eq. (|6| and 
d7l), are practically the same (see Figure [3]), except when fi < lMeV. In 
this case, the last two terms of Eq. ^ arising from the q^q u / part of the 
propagator of the axial vector particle start to dominate the amplitude. 

Our bounds on the parameters g v and g\ (Figure [3| are rather strong; 
say, for \i = 10 MeV, g VA < 5 ■ 10~ 3 at 90% C.L., which corresponds to 

g V N A < 0-071 , (8) 

four times smaller than the QED coupling constant y/A-na ~ 0.3. For the 
scalar exchange, taking fi = 10 MeV, we get a much weaker bound, g 2 s < 1.4. 



3 Bounds from lower energy data 

We will now compare our results from the previous section with other searches 
for new interactions [5] [H] in which new light particles participate. In the 
literature, the effect of new forces is usually parameterized as a deviation 
from the Newtonian gravitational potential: 

V (r) = -G N r ^[l + a G exp(-r/X)] , (9) 

which is an adequate approximation for the description of the effect of a new 
particle exchange between nonrelativistic constituents. The following rela- 
tionship exists between the coupling constant «g and characteristic length A 
and our parameters gf and fi in cases of vector and scalar exchanges: 

Q 2 

a ° = = 135 ' 1037 ^ • lgo = lg9 ^ + 3713 ■ < 10 > 

A(<m) %(NfeV)LE.10" - ■gA=-lg/.-10.7 ■ (ID 

The pseudoscalar exchange would not modify the potential in a nonrel- 
ativistic approximation, while axial coupling leads to an interaction among 
spins of constituents. 

(A) References [5]-[7] show e.g. that values of «g larger than 1 are 
excluded for A larger than 0.1 mm, while for smaller A the upper bounds 
on «g rapidly grow, reaching 10 9 at the micron scale. 
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In papers [E] [H], analogous bounds for shorter distances are presented. 
We see e.g. that for A = 1CT 13 m, «g should be less than 10 30 , or g vs less 
than 10~ 7 . The corresponding value of \i is 2 MeV. In Figure |4j our limits 
for the V and S particle exchanges in terms of and A parameters are 
compared with the limits obtained in papers |5]-|14j. 

(B) Additionally, the data on low-energy (1 keV < E n < 10 keV) neutron 
scattering on 208 Pb [15] were applied in paper [16] to obtain bounds on the 
possible contributions of a light scalar particle exchange to neutron-nucleus 
potential. The upper bound on the coupling constant of a 10 MeV boson to 
a nucleon obtained in [TS] corresponds to g\r S < 4 • 10~ 6 . This rather restric- 
tive bound was obtained from the analysis of the shape of the differential 
cross-section of low-energy n — 208 Pb scattering, where the additional term, 
originating from the light scalar boson exchange, leads to a modification of 
angular dependence not observed in the experimental data. Let us stress 
that the same bound is valid for g v . 

According to [13 US], experimental data in the keV energy range are well 
described by the following expression: 

^ = ^[l + ^co S 9] , (12) 

where ^(Jq/Att ps 10 fm, and uj = (1.91 ± 0.42)10~ 3 keV -1 . These numerical 
values are very reasonable from the nuclear scattering point of view and, 
from the demand that these values are not spoiled by a Yukawa potential 
contribution originating from a light boson exchange, bounds on g and /i were 
obtained in [16]. The point is that the Yukawa amplitude, interfering with 
the strong interaction amplitude, will show up in the following contribution 
to uj for E — > 0: 

|Ao,| = J^L^4 , (13) 



and from the demand that Auj < uj, the above mentioned bound was ex- 
tracted. For an update of results obtained in [16], see [PT] . 

(C) It is quite natural to assume that the coupling of a new light boson 
with nucleons originates from its coupling with u- and (i-quarks. In this case, 
bounds from pion and kaon decays [TS] are applicable. Let us start with 
vector coupling. According to CVC, couplings to nucleons are equal to the 
sum of the couplings to quarks: 2f u y + fdv for a proton and f u y + Ifay for 
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Figure 4: Experimental limits on olq and A from [5]-[T3] parameterizing devia- 
tions from Newton's law. Our limits transformed into coordinates olq and A are 
also shown for comparison. 
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a neutron (/$ are analogous to our g l N ). The 7r° — > VV decay contributes to 
7T° — > invisible decays and, using the experimental bound Br(7t° — > i/v) < 
2.7 ■ 1CT 7 in [18J, the following bound was obtained: 

\l\flv~fh\ < 4 ■ > (14) 

which is automatically satisfied for an isoscalar coupling, f u y = fdv- How- 
ever, the bound on the tt° — > 7V decay, which contributes into 7r° — > 
7 + invisible mode, allows bounding of isoscalar couplings as well [18] : 



2fuV + fdV <16 1Q -3 (15) 

Here, the experimental bound Br (it — > 71/1/) < 6 • 10~ 4 was used. These 
numbers should be compared with our bound on (III). 

Since 7r° — > S'S and 7r° — > 5*7 decays violate the corresponding P- and 
C-parities, we do not obtain bounds on fs from these decays. C-parity 
conservation forbids the tt° — ► 7A decay as well, while from the bound on 



7T° — > invisible decays, we get the coupling constant bound ( 14 ) for the axial 
vector boson. 

More stringent upper bounds on the coupling constants follow from very 
strong experimental limits on the branching ratio Br(K + — > 7r + + vv) < 
2 • 1CT 10 . The longitudinal component of the axial vector boson contributes 
to the decay amplitude proportionally as (2m q //j,)f q A [IE], and even if the 
axial vector boson couples only with light quarks, we obtain: 

fu,dA < KrV(MeV) . (16) 

The factor of 2m q / fi is absent when the axial vector interaction is substi- 
tuted by the scalar interaction, and thus we obtain: 

fu, d s < 10- 5 . (17) 

Fortunately, CVC forbids K — >• decays for yU 2 = 0, so that is why the 
bound on the vector coupling for light \x is not very strong: 

fu, d v(— V <10~ 5 • (18) 
\m K J 
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4 Conclusions 



Our bounds obtained from high-energy and very small momentum transfer 
np elastic scattering data [2] provide exclusions of new forces at distances 
above 5 fermi, which corresponds to exchanged particle masses lighter than 40 
MeV. These bounds are extracted in a covariant approach, as an alternative 
to the bounds on couplings at larger distances, extracted from the absence 
of deviations from the Newtonian gravitational law. 

Both low-energy n — 208 Pb and high-energy np scattering data lead to 
similar upper bounds on the coupling constants for w 10 MeV vector bosons, 
though upper bounds from n - Pb scattering on the coupling constant are 
~30 times lower and close to the bounds from 7r° — > invisible and 7r° — > 7 + 
invisible decays on the vector coupling constants with quarks. 

Strong upper bounds on the "new physics" contribution into the K + — > 
7r + + invisible decay allows us to get very strong bounds for scalar and axial 

AS ^ 

vector bosons: g N ' < 10 for a 10 MeV boson mass. 

We are grateful to A. B. Kaidalov, B. Z. Kopeliovich, L. B. Okun, and Yu. 
N. Pokotilovsky for useful discussions and comments. M. V. was partially 
supported by Rosatom and grants RFBR 07-02-00021, RFBR 08-02-00494 
and NSh-4568.2008.2. Y.K. would like to acknowledge research support from 
the United States Department of Energy HEP grant DE-FG02-91ER40627. 



5 Appendix 

For vanishing light quark masses, their isotriplet axial current is conserved. 
The same should hold for the nucleon currents and is achieved by accounting 
for pion exchanges: 



A, 



k a k 



'ft 



k 2 — m 2 



(9, 



a/j 



k 2 — ji 2 



9p.u 



knky 



mi 



9l 



k 2 -u 2 



k a kf3 {m\ - 2jjL 2 m 2 T + k 2 fi 2 ^ 



9a(3 



(k 2 — m 2 ) 2 



k 2 



(19) 



For a massless pion, the 1/fJ 2 singularity cancels out and the expression 
in square brackets contains k a k^/k 2 , which, acting upon fermionic axial cur- 
rents, becomes (2mjv) 2 /A; 2 . The numerator of the expression for differential 
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cross-section is regular at k 2 = t = since the square of the pseudoscalar 
exchange amplitude contains t 2 in the numerator (see (5)), while the inter- 
ference of the axial vector and pseudoscalar exchanges is proportional to t 
(the denominator equals (t — fi 2 ) 2 independently of the spin of the exchanged 
boson). 

In real life, light quarks, as well as pions, have nonzero masses, and to 
obtain an amplitude square for the axial vector boson exchange we should 
substitute \i by \x in the square brackets of (7), where 

1 1 ml — 2fi 2 ml + tfi 2 
jl 2 fi 2 (t — m 2 ) 2 

and for t, fi 2 <C m 2 we get jl — /i. 

The numerator of the expression for the differential cross-section is sin- 
gular for fi — > 0. However, in renormalizable theory, the mass of the axial 
vector boson equals its gauge coupling constant (g^ in our case) times the 
vacuum average of the corresponding higgs field. 

As an example, one can have in mind the expansion of the Standard Model 
with two higgs doublets with opposite hypercharges, where Peccei-Quinn 
[/(l)-symmetry is spontaneously broken producing an axion. In order to 
suppress axion couplings to quarks and leptons, the additional singlet neutral 
higgs field N is usually added, which makes the axion invisible. Gauging of 
Peccei-Quinn U(l) leads to the axial coupling of the corresponding vector 
boson to matter. Such a light axial vector boson is discussed in particular 
in [18] . where it is light due to the smallness of the gauge coupling constant, 
while the vacuum average < N > ^> 100 GeV, making it superweakly 
coupled to matter (gA./fJ> ~ 1/ < iV >). 
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